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ABSTRACT 



Context. S147 is a large faint shell-type supernova remnant (SNR) known for its remarkable spectral break at cm- 
wavelength, which is an important physical property to characterize SNR evolution. However, the spectral break is 
based on radio observations with limited precision. 

Aims. New sensitive observations at high frequencies are required for a detailed study of the spectral properties and the 
magnetic field structure of S147. 

Methods. We conducted new radio continuum and polarization observations of S147 at All cm and at A6 cm with the 
Effelsberg 100-m telescope and the Urumqi 25-m telescope, respectively. We combined these new data with published 
lower-frequency data from the Effelsberg 100-m telescope, and with very high-frequency data from WMAP to investigate 
the spectral turnover and polarization properties of S147. 

Results. S147 consists of numerous filaments embedded in diffuse emission. We found that the integrated fiux densities 
of S147 are 34.8 ± 4.0 Jy at All cm and 15.4 ± 3.0 Jy at A6 cm. These new measurements confirm the known spectral 
turnover at ~1.5 GHz, which can be entirely attributed to the difi'use emission component. The spectral index above 
the turnover is a = —1.35 ± 0.20 (S ~ 1/°'). The filamentary emission component has a constant spectral index over the 
entire wavelength range up to 40.7 GHz of a = —0.35 ± 0.15. The weak polarized emission of S147 is at the same level 
as the ambient diffuse Galactic polarization. The rotation measure of the eastern filamentary shell is about —70 rad 
m~^. 

Conclusions. The filamentary and diffuse emission components of S147 have different physical properties, which make 
S147 outstanding among shell type SNRs. We attribute the weak polarization of S147 at All cm and at A6 cm to a 
section of the S147 shell showing a tangetial magnetic field direction. 

Key words. ISM: supernova remnants - radio continuum: ISM - techniques: polarimetric 



1. Introduction 

The extended object Shajn 147 (S147) has long, beauti- 
ful delicate filaments visible in the optical bands, and is 
located in the anti-center of the Galaxy. It was first classi- 
fied as a p o ssible shel l-type supernova remn ant (SNR) by 
iMinkowskil (|1958D and Ivan den Berghl (Il960l). The distance 
to S1 4 7 was estimated to b e around 1 kpc ("ciark fc Caswelll 
119761: iKundu et al.| [T980l). using the so called E-D rela- 
tion (iMilne I 1979 ). The interstellar reddening toward the 
SNR ( Fesen et al.l [TMl) suggests a smaller distance of 
0.8 kpc. This is supported by absorption hues of the Ble 
star HD3666 5 at 880 pc originating from high velocit y 
gas of S147 (|PhilliDS et al.l 119811 : ISallmen fc WelshI l2004f) . 
The expansion velocity of the shell is about 80 - 120 km 
s-i (jPhillips et aTlll98ll : lKirshner fc ArnoTdlll979D . The age 
of S147 was estimated to be 2 x 10^ yr through the 
radiative blast wave's S edov- Taylor solution ( Sofue et al.l 
119801: IKundu et al.lll980D . The PSR J0538 +2817 has been 
disco vered within the boundary of S147 (| Anderson et al.l 
Il996l ). and is probably physically associated with the SNR 
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* Figs. 1 and 2 are available in FITS format at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A-|-A/| 



()Anderson et al.lll996t iNg et al.ll2007D . The kinematic age 
of the pulsar, if originated from the SNR center, is about 
^ 4 X 10^ yr. Th e parallax of the pulsar suggests a distance 
of lA7t°o f^ kpc (INg et al.ll2007l) . 

iDeNoveil (|l974l ) summarized early radio maps of S147 
between 178 MHz and 1420 MHz. She obtained a flux den- 
sity spectral index of a = —0.67 ± 0.37, revealing S147 
as a SNR. Later, radio obser vations of S147 at 1 648 MHz 
and 2700 MHz were made bv lKundu et all |1980), and the 
southern part of the SNR was observed at 4995 MHz by 
ISofue et all (|1980[ ). These high-frequency observations sug- 
gested a spectral turnover near 1.5 GHz with a flat spec- 
trum at lower frequencies and a steep spectrum at higher 
frequencies. The turnover needs to be confirmed due to 
large uncertainties in the low-frequency measurements and 
the fact that the observation at 4995 MHz covered only 
the southern part of the remnant. Later sensitive mea- 
surements, at 1425 MHz and 2695 MHz, by Fiirst fc ReicB 
(|1986D could not improve the integrated radio spectrum, 
but revealed spectral index variations across S147. Regions 
of bright filaments show a flux density spectral index of 
about a = —0.5, while regions in between show a spec- 
tra l index near a = — 1.0. The observation at 863 MHz 
by iReich et all (|2003t) provided a more precise flux den- 
sity at a low frequency, but did not improve on the un- 
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Table 1. Observational parameters 



Wavelength 


A6 cm 


All cm 


Frequency 


4800 MHz 


2639 MHz 


Bandwidth 


600 MHz 


80 MHz 


HPBW ['] 


9.5 


4.4 


aperture efficiency [%] 


62 


53 


beam efficiency [%] 


67 


58 


Tsys[K] 


22 


17 


Ts[i^]/S[Jy] 


0.164 


2.60 


Main Calibrator 


3C286 


3C286 


Flux Density 


7.5 Jy 


11.5 Jy 


Polarization Percentage 


11.3% 


9.9% 


Polarization Angle 


33° 


33° 


No. of coverages 


14 


6 


pixel integration time [sec] 


18.6 


3 


r.m.s. (total intensity) [mK] 


0.7' 


5.0 


r.m.s. (polarized intensity) [mK] 


0.3* 


3.5 



*The quoted sensitivity at A6 cm 
1=178?0 to 1=182?2, b=-3?6 to b= 



is valid for the area 
=-0'?6 (see text) 



cer tainty in the spectrum . The 2695 MHz radio map made 
bv lFiirst fc ReichI (jl986f ) for the first time showed Unearly 
polarized emission associated with some filamentary struc- 
ture. Detailed information from this polarization observa- 
tion cannot be derived due to an insufficient signal-to-noise 
ratio. 

One goal of this work is to use the new sensitive A6 
cm receiver at the Urumqi 25-m radio telescope to obtain 
a complete radio map at this band for a detailed study 
of the spectral turnover. Here we present results of to- 
tal power and linear polarization maps at 4800 MHz (A6 
cm). We also present new observations at 2639 MHz (All 
cm) with the Effelsberg 100-m radio telescope and unpub- 
lished polarization data at 1420 MHz (A21 cm) obtained 
from the "Effelsb erg Medium Galactic Latitude Survey" 
(jReich et al.ll2004[ ). The new observations and the data pro- 
cessing are described in Sect. 2. The investigation of the 
integrated radio spectrum and possible spectral variations 
across the remnant are discussed in Sect. 3. The radio polar- 
ization and the rotation measure toward S147 are analyzed 
in Sect. 4. The results are summarized in Sect. 5. 



2. Observations and data processing 

New observations of S147 have been made with the Urumqi 
25-m radio telescope at A6 cm and Effelsberg 100-m radio 
telescope at All cm. The main parameters are listed in 
Table [II 



2.1. Observations at X6 cm 

We made continuum and polarization observations of S147 
between January 2005 and February 2006. We scanned the 
large area of S147 along Galactic longitude or Galactic lat- 
itude direction. Five coverages of the field of {Al x Ab) = 
(5° X 5°) centered on {I, h) = 180?2, -1?7 were observed in 
January 2005, four coverages (A^ x Ab) = (4?2 x 4?2) cen- 
tered on {I, b) = 180?1, -1?5 were added at the end of 2005, 
and five additional maps of {Al x Ah) — (5° x 5°) centered on 
{I, b) = 180?2, — 1?6 were observed again in February 2006. 
We moved the telescope with a scanning velocity between 
2'/s and 3'/s. The sampling and the pixel size of the map 
were 3'. There is one map scanned with a narrow bandwidth 



(295 MHz instead of 600 MHz) to mask interference from 
InSat (see lSun et al.l[2007l ). so the effective total integration 
time is about 18.6 s per pixel and the theoretical sensitivity 
is about 0.2 mK Ta, using ct/ — Tgys/V Ai/t, where Tgys is 
the system temperature, Ta is the antenna temperature, 
Az^ is the bandwidth, and r is the integration time. The 
quoted sensitivity corresponds to 0.3 mK Tb, where Tb is 
the brightness temperature. 

The raw data from the four backend channels (RR*, 
LL*, RL*, LR*), together with time information and 
telescope position were stored in the MBFITS format 
(jMuders et al.ll2007D. and sub sequently converted into the 
N0D2 format ( HaslamI |1974| ) to be further processed by 
the N0D2 based data reduction package under a Linux 
operation system. We edited the individual maps to re- 
move any interference, to correct baseline curvatures by 
polynomial fitting and to suppress "scann ing effects" by 
applying an "unsharp masking method" (|Sofue &: ReichI 
|l979) especially for the U and Q maps. Finally, we re- 
tabulated all maps into a map size of 5° x 5° centered on 
l,b= 180?2, -1?7 with a pixel size of 4' and "weaved" all 
maps together applying the method of lEmerson fc Gravi 
()l988l ). which "destripes" a set of maps in the Fourier do- 
main. We finally corrected the polarization map from t he 
positive noise bias following I Wardle fc Kronberg I (|l974f ). 

The Urumqi A6 cm total intensity map of S147 and po- 
larization intensity map with vectors overlaid in magnetic 
field direction are shown in Fig. [TJ The total intensity map 
clearly reveals the spherical SNR shell with a hollow re- 
gion in its center. The r.m.s. -noise in the common area of 
all coverages (see Table [1]) in total intensity is 0.7 mK Tb 
and in polarized intensity is 0.3 mK Tb. These values are 
slightly la rger than expect e d. It c ould be partially caused by 
confusion. lUvaniker et al.l (|l99! /l found a confusion limit of 
15 mK Tb at 1400 MHz, which converts to 0.48 mK Tb at 
A6 cm based on a flux density spectral index of a = —0.9 
(jZhang et al.ll20d3[ ). Combining this value with the theo- 
retically expected value, we obtain an r.m.s. -noise of 0.5 
mK Tb. Another part may be explained by residual gain 
fluctuations, "scanning effects" , atmospheric and ground 
radiation v ariations, and lo w-level interference close to the 
noise level (|Sun et al.l 120061 ). 



2.2. Observations at Xll cm 

We observed S147 at All cm with a new receiver in- 
stalled in 2005 in the secondary focus of the Effelsberg 
100-m radio telescope. The map size is 5° x 5° centered at 
l,b = 180?2,— 1?6. We made six coverages along Galactic 
longitude and Galactic latitude direction with a scanning 
velocity of 4'/s. We processed the maps in the N0D2 format 
based standard reduction package for observations with the 
Effelsberg 100-m telescope. Bad data caused by strong in- 
terference were blanked. Scanning effects were suppressed 
by using the method of "unsharp masking" developed by 
ISofue fc ReichI (|l979t) . Finally, we "wea ved" together all 
edited maps applying the method of E merson fc Gravj 
(Il98l . 

The All cm maps are shown in Fig. [51 The r.m.s.- 
noise in total intensity (upper panel) is measured to be 
5.0 mK Tb. In polarized intensity, we found a value of 3.5 
mK Tb. Because of the higher angular resolution of 4f4, 
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Table 2. Strong point-like sources in the boundary of S147 detected by Effelsberg at All cm and Urumqi at A6 cm. The 
flux density is given by the peak flux of two-dimensional gaussian fits. 



Source 


ai950 


5l950 12000 


(52000 








(h m s) 


(° ' ") (h m s) 


(° ' ") 


(mJy) 


(mJy) 



0531+275 


05 31 18.3 


27 30 24 


05 34 26.6 


27 32 23 


533+37 


287+44 


0531+279 


05 31 21.7 


27 53 41 


05 34 30.6 


27 55 39 


383+33 


265+24 


0533+282 


05 32 48.2 


28 09 24 


05 35 57.5 


28 11 16 


135+11 


51+14 


0533+272 


05 33 05.6 


27 10 43 


05 36 13.5 


27 12 33 


115+13 


84+13 


0534+284 


05 34 13.5 


28 24 41 


05 37 23.2 


28 26 26 


107+13 


130+20 


0535+266 


05 34 33.1 


26 37 25 


05 37 40.2 


26 39 09 


86+12 




0536+285 


05 36 20.5 


28 30 19 


05 39 30.4 


28 31 56 


133+10 


71+18 


0538+287 


05 38 06.3 


28 41 14 


05 41 16.5 


28 42 43 


762+46 


527+47 


0539+290 


05 39 02.7 


29 00 20 


05 42 13.4 


29 01 45 


665+42 


435+33 


0541+269 


05 41 03.4 


26 55 07 


05 44 11.0 


26 56 23 


146+16 




0544+273 


05 44 27.5 


27 20 58 


05 47 35.7 


27 21 59 


262+17 


251+22 





Fig. 3. Left panel: The filamentary structure of S147 at All cm with most of strong point-like sources and diffuse 
large-scale (> 12') emission component subtracted (see text for details). Right panel: the optical image of S147 (Credit: 
Digitized Sky Survey, ESA/ESQ/NASA, http://antwrp.gsfc.nasa.gov/apod/ap051129.html). 



the shell and the filamentary structures can be more clearly 
identified than in the A6 cm map. 



3. Comparison of radio and optical observations 

We used All cm total intensity map at the original angular 
resolution after subtraction of the point-like sources listed 
in Table [5] to separate small-scale and large-scale structures 
applying the "background filter ing" (BGF) procedure, in- 
vented bv lSofue fc ReichI ()1979[ ). where a filtering beam of 
12' was used. A comparison of the small-scale radio struc- 
ture with optical filaments of S147 is displayed in Fig. [31 
which demonstrates an excellent agreement, even faint opti- 
cal structures of the blow-outs toward north-east and south- 
west are visible in both bands. 

We also overlay the radio filamentary structures on the 
Ha map compiled by Finkbeiner (2003) in Fig. IH and find 
a rather good positional coincidence when inspecting in- 
dividual filaments. However, we find large variations in 
th e relative int e nsitie s, which have been already mentioned 
bv ISofue et al.l ()198G[ ). These variations may be explained 
by either stronger magnetic fields in the radio bright fila- 
ments and/or lower temperatures/densities in the optically- 
bright filaments. The measured temperature and density 



variations (jFesen et al.lll98 5': D'O dorico &: S abbadin' 19771) 
support this conclusionT lP'Odorico fc Sabbadin. ((1977. ) ob- 
tained an average electron density of about 250 cm~^ as- 
suming an electron temperature of 10^ K and showed that 
the compression by an adiabatic shock is insufficient to ob- 
tain these densities. Obviously, at least parts of S147 are 
already in the cooling phase of SNR evolution. 

4. Spectral properties of S147 

4.1. The integrated radio spectrum of S147 

The new observations with the Urumqi 25-m radio tele- 
scope provide the first complete map of S147 at A6 cm. In 
order to estimate the integrated flux density of S147, we 
subtracted eleven strong point-like sources (see Table [5]) 
toward S147 from the A6 cm and All cm maps. The flux 
densities of these sources have previously been listed by 
iFiirst et al.l (|1982f ). and agree within 0.1 Jy to those in 
Table [H We calculated the integrated flux density with the 
method of ring integration by summing up the emission in 
concentric rings centered on l,b = 180?2,— 1?6. The inte- 
gration stopped at a radius of 2?0. Beyond this radius, we 
found an average base-level of 8.7 mK Tb at All cm and 
of 1.9 mk Tb at A6 cm. We have subtracted these base- 
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Fig. 1. The Uriimqi A6 cm map of S147 at an angular 
resolution of 9.'5. The upper panel shows the total inten- 
sity map. The lower panel shows the polarization intensity 
map with contours of total intensities overlaid (the strong 
point-like sources listed in Table 2 have been subtracted). 
Contours start at 10 mK Tb, and increase by a factor of 
1/2. The bars show the orientation of the magnetic field 
(E-|-90°) for the case of negligible Faraday rotation. Their 
length is proportional to the polarized intensity. 



levels from the maps before calculating the flux density of 
S147. Taking these local base-levels into account, polygon 
integrations at the two wavelengths just outside the pe- 
riphery of SI 47 result in the same flux density values as 
obtained from the method of ring integration. From vari- 
ations just outside S147 we found the uncertainty in the 
base-levels to be 2 mk Tb at All cm and 1 mk Tb at A6 
cm. Considering 5% of other uncertainties including cali- 
bration, the flux density values are 34.8 ± 4.0 Jy at All 
cm and 15.4 ± 3.0 Jy at A6 cm. Figure [H] displays the inte- 
grated spectrum of S147 including available low-frequency 
flux density data. 



mi 
-a 

?. -r 



« 

• 

• 




m 

* 

• 

■ 


* 


■ 




* 






■ 

■ 


• 


* 

^ « 



'62' ler 'sc- i73« 
Galaciic Lonq'tjde 
Dfelsberq ^' 1 cn P -t B 



176" 



3 -r 




132" -61° ISO' 179' 

Galactic Lsngituda 



Fig. 2. Same as Fig. [l]but showing the Effelsberg All cm 
maps at an angular resolution of A'A. Contours start from 
50 mK Tb and increase by a factor of 2. 



iReich et al.l (|2003[ ) placed the turnover frequency be- 
tween 1.5 GHz and 2 GHz using the extrapolat ed flux den- 
sity of the incomplete 5 GHz map observed by ISofue et al.l 
(|1980D . The extrapolated flux density value as well as the 
spectral turnover at ^ 1.5 GHz are clearly confirmed by 
new measurements at A6 cm and All cm. The spectrum 
is fiat at lower frequencies and steepens toward frequen- 
cies higher than 1 GHz to 2 GHz. The low-frequency spec- 
trum has a spectral index of about a ^ —0.30 ± 0.15. The 
spectral index is about a ~ — 1.20± 0.30 above the break 
frequency. The relatively large error at high frequencies is 
mainly caused by the still large uncertainty of the flux den- 
sity at 4800 MHz (A6 cm), which, in principle, is very dif- 
ficult to reduce for faint large diameter objects like S147 
located in the Galactic plane. 
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Fig. 4. The filamentary str uctures of S147 ob served at All 
cm is overlaid with Hq, map (|Finkbeinerj|2003l) . The angular 
resolution of the Hq map is 6'. The contours at All cm start 
at 20 mK Tb and run in steps of 35 mK Tb- 
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Fig. 5. Spectrum of the integrated radio flux densities of 
S147. The low-frequency spectrum has a spectral index of 
a ^ —0.30 ± 0.15. Above the break the spectral index is 
a ~ —1.20 ± 0.30. The S147 inte grated flux d ensities are 
taken from lFurst fc ReichI (|1986D ; iReich et all (|2003D ; and 



iKovalenko et al.1 ( 1994D. We q r nitted the A6 cm flux density 
value obtained bv lSofue et al.l (|l980f) . be cause this in e asure - 
ment did not cover the entire source. iReich et al.l (|2003D 
have corrected the values for the contr ibution of point-like 
sources as listed by [Fiirst et al.l (|1982D . assuming a mean 
source spectral index of a = —0.75 at frequencies lower 
than 408 MHz. 



4.2. The spectral index map of S147 

The Urumqi A6 cm and the Effelsberg All cm maps with 
point-like sources and constant base-levels removed were 
both convolved to a common angular resolution of 10'. Wc 
estimated the relative position accuracy of the two maps 
by using the positions of point-like sources from Table 2. 
We found the mean positional difference to be 6" with a 
maximum of 28". 
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Fig. 6. Spectral index map calculated between All cm and 
A6 cm at 10' angular resolution. The overlaid contours rep- 
resent total intensities at All cm (with point-like sources 
listed in Table 2 subtracted), starting from 50 mK Tb and 
increasing by a factor of 2. 



We display the spectral index map between All cm and 
A6 cm in Fig. [SI We calculated the spectral index of each 
pixel from the brightness temperatures at the two frequen- 
cies. In order to achieve accurate spectral indices and to 
exclude the influence of systematic effects we set a lower 
flux density limit of 40 mK Tb and 6 mK Tb for the All 
cm and A6 cm map, respectively. The possible variations of 
the base- levels at All cm and A6 cm cause an uncertainty of 
the spectral indices of Aa ~ 0.3. The uncertainty is largest 
where the total intensity is small. The general result of a 
much flatter radio spectrum of the filaments compared with 
the large-scale component is not altered. However, system- 
atic errors make it difficult to derive positional spectral 
index variations of the large scale component. The spec- 
trum for regions with bright radio filaments, particularly 
in the eastern and northern part of S147, is much ffatter 
than for regions with less pronounced filaments and for the 
large-scale diffuse emission component near the center. As 
seen in the total intensity maps, faint, extended Galactic 
emission exists outside of S147. Such emission features may 
also exist in the projected region toward S147, leading to 
variations of the spectral index. 

4.3. T-T plot analysis 

The temperature— temperature (T-T) plot is a method 
([Turtle et al.lll96^ to investigate the spectral behavior, in- 
dependent from a constant base-level across the source. We 
used the maps at an angular resolution of 10' at All cm 
and A6 cm to check the spectral index of the filamentary 
and diffuse emission components. We first decomposed both 
maps into small-scale filamentary structures and large-scale 
diffuse emi ssion component with the BGF procedure de- 
scribed by ISofue fc ReichI ([1979f ) with a filtering beam of 
18' and cut-values of 100 mK Tb at All cm and 25 mK Tb 
at A6 cm. The resultant TT-plots for the filaments and for 
the diffuse emission component are shown in Fig. [7| They 
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Fig. 7. The T-T plot for the filamentary and diffuse emission components of S147 between All cm and A6 cm. Filamentary 
emission component in general has a spectral index of a = —0.43 ± 0.10. The diffuse emission component has a spectral 
index of a = -1.35 ±0.20. 



have largely different spectral indices. We found the steep 
spectrum with a spectral index of a = —1.35 ± 0.20 for 
the diffuse emission component, while we found a spectral 
index a = —0.43 ± 0.10 for the filamentary emission com- 
ponent. This index is slightly steeper than the integrated 
low-frequency spectrum (a — — 0.30±0.15), but is compara- 
ble considering the quoted uncertainties. The quoted errors 
include the statistcal error (number of data points reduced 
to one per HPBW) as well as the error, which follows from 
fitting the data twice, alternatively taking the data of one 
of the two wavelengths as independent variable. The result 
in Fig. [7] may be affected by possible variations of the base- 
levels across S147. This was inspected using TT-plots of 
four smaller regions separated by 1=180?0 and b=-l?9. For 
the filamentary component, the variation of the spectral in- 
dex is within the quoted error of Aa = ±0.10. Therefore, 
the filamentary emission component probably has a fairly 
straight radio spectrum over the accessible frequency range. 
For the diffuse emission component the spectral index varies 
between q:=-1.10 and q;=-1.56. There is no doubt that the 
radio spectrum of the diffuse emission component is much 
steeper than that of the filamentary component. This is in 
agreement with the result from the spectral index map. It 
is entirely the diffuse emission component that produces 
the spectral turnover around 1.5 GHz. The integrated flux 
densities of the separated diffuse and filamentary emission 
components are ~ 22.8 Jy and ^ 12 Jy at All cm, and 
^ 8.4 Jy and ~ 7 Jy at A6 cm, respectively, which also 
shows the difference of spectral indices of the two compo- 
nents. 

It is clear that above 4800 MHz the flux density of the 
radio filaments dominates the integrated flux density of 
S147, causing a spectral flattening toward higher frequen- 
cies. To check this, we inspected the WM AP data in all 
five bands between 22.8 GHz and 93.5 GHz (Hinsha w et al.l 
[2007 ) for a signature of S147. We are able to clearly trace 
S147 up to the highest frequency especially its region with 
strong radio filaments near I = 181°. 
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Fig. 8. The filamentary emission component of S147 at 40.7 
GHz (grayscale) observed by WMAP at an angular resolu- 
tion of 49'. All cm contours of Effelsberg measurements at 
an angular resolution of A'A start at 20 mK Tb and run in 
steps of 35 mK Tb. 



We convolved all WMAP data and the All cm map to 
the beam size of the 22.8 GHz map (49'). We separated 
the filamentary emission component of SI 47 in all maps by 
applying the BGF procedure described by ISofue &: ReichI 
(1979), using a filtering beam size of 50'. As an example, 
we show the WMAP data at 40.7 GHz in Fig.[5]at 49' angu- 
lar resolution in comparison with the filamentary emission 
component of the Effelsberg All cm map at the original an- 
gular resolution. To check the spectral index, we performed 
a T-T Plot (Fig. [5]) for the filamentary emission component 
at 40.7 GHz and All cm of the area I = 180?8 to / = 181?9, 
b = -0?87 to 6 = -2?7. We found the spectral index to be 
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Fig. 9. The T-T Plot of filamentary emission component 
at A7 mm (40.7 GHz) and All cm. The spectral index is 
a = -0.31 ±0.20. 
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Fig. 10. Integrated total intensity spectrum of S147, in- 
cluding WMAP data for the filaments. The weaker con- 
tribution from the diffuse emission component can not be 
obtained from the WMAP data. The model spectrum of 
the filamentary component, the diffuse component with a 
spectral break at 1.7 GHz, and the sum of both are shown 
in addition (see text). 



a = — 0.31±0.20, which is consistent with that between All 
cm and A6 cm. A similar spectral index for the filamentary 
emission component is obtained between All cm and 22.8 
GHz and 33.0 GHz, respectively. 

From the various methods we may derive a spectral 
index for the filamentary emission component of a = 
—0.35 ± 0.15. Using the flux density of 7 Jy for the fila- 
mentary structures at A6 cm, we calculate the flux density 
of filamentary emission component to 3.5 Jy at 40.7 GHz, 
4.3 Jy at 22.8 GHz, and 3.8 Jy at 33.0 GHz. The fiux den- 
sity of the diffuse emission component of S147 at 40.7 GHz 
or other WMAP frequencies can not be determined since a 
separation of this component from the underlying stronger 
large-scale Galactic emission is not possible. 

At 60.8 GHz and 93.5 GHz the emission in the quoted 
area increases significantly, which is most likely due to ther- 
mal dust emission associated with the SNR. 



We show all fiux densities up to 40.7 GHz in Fig. (TUl 
where the WMAP values miss the small contribution from 
the diffuse emission component. In addition, we have plot- 
ted the supposed straight spectrum of the filamentary com- 
ponent with a spectral index of a — —0.35. We have also 
plotted the contribution of the diffuse large-scale emission 
component and the sum of both. The turnover of the inte- 
grated total intensity spectrum has to be attributed to the 
turnover of the large-scale emission component. We assume 
that both emission components have the same spectral in- 
dex below the break frequency. In summary, the integrated 
total intensity spectrum is flat at low frequencies, steep be- 
tween 1 GHz and about 10 GHz, and flattens again toward 
higher frequencies. No other SNR is known showing a sim- 
ilar radio spectrum. Hence, the radio spectrum of S147 is 
unique. 

4.4. The analysis of the spectral turnover 

A spectral turnover is very rare for shell-type SNRs. Some 
early claims have turned out to result from limited accuracy 
of high-frequency measurements. S147 is the best known 
case for a spectral break so far. 

Several possible mechanisms cause a spectral turnover: 
(1) The presence of two populations with different energy 
spectra; (2) shock acceleration of particles; (3) the com- 
pression of the Galactic magnetic field; and (4) synchrotron 
losses (aging) of electrons^ 

(1), (2) According to lLeahv fc Rogej (|l998l ) the first two 
mechanisms may not be relevant for a synchrotron spec- 
trum as given in Fig. [51 

(3) The Galactic relativistic electron spectrum has a 
slow turnover due to various loss and propagation effects 
over time. .WebsteS (1974) has compiled the Galactic radio 
spectrum, which varies from a ~-0.5 to a ~-0.8 between 
'^lOO MHz and ~600 MHz. The corresponding electron en- 
ergy index varies from ^2.0 to ^2.6 over the energy range 
2.6 GeV to 6.5 GeV in case the magnetic field strength is 
3 /iG. Direct observations of the electron spectrum ( Tan3 
11984ft also reveal a change of the Galactic electron spectral 
index from 2.0 at low energy (< 2 GeV) to 3.5 at higher en- 
ergies. Due to the shock wave of the SNR, the compressed 
local magnetic field in a SNR shifts the spectral turnover 
towar d a highe r frequ ency, e.g., that at about 1.5 GHz of 
S147. iDeNoveij ()1974D has studied this case extensively. A 
crucial parameter is the radio flux density, which is re- 
lated to the observed turnover frequency and the distance. 
Following the argument by DeNoyer, it is possible to ex- 
plain the observed radio spectrum of SI 47, in particular, 
for the smaller distance of <1 kpc. 

(4) If synchrotron losses of electrons occur uniformly 
over the whole lifetime of the SNR, the turnover would 
happen at :/c(MHz) = 3.4 x 10^ B(10//G) t(lOV) 
For S147 the turnover frequency i^c is observed to be at 
1.5 GHz and its age is about lO^yr, so the magnetic fields 
should be about 0.3 mG. Such a high magnetic field is not 
known to exist in shell type SNRs. DeNoyer (1974) dis- 
cusses the situation that the acceleration and synchrotron 
losses occur during the early phase of SNRs. iGullI ()1973l ) 
has shown that in the early phase Rayleigh- Taylor instabil- 
ities create a convective zone, which can sustain magnetic 
fields B of 1000 /iG within young remnants for approxi- 
mately 100 years, corresponding to a radius Rc of the rem- 
nant of 1 pc in the low-density region of the Galactic 
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anti-center. Due to the subsequent expansion, the mag- 
netic field decreases with (R/Rc)~^, while the energy de- 
creases with (R/Rc)~^. The turnover frequency decreases 
with (R/Rc)-^. For S147, R/Rc is 21 for a distance of 800 
pc corresponding to a turnover fr equency of 12 M Hz and 
154 MHz, respectively. However, as lDeNoveil ()1974| ) pointed 
out, at this stage of evolution the magnetic field of the SNR 
is dominated by the shock compressed interstellar magnetic 
field, which may be about 10 /xG in the large-scale emis- 
sion area. In that case, the turnover frequency decreases 
with (R/Rc)-2 X (10 /xG/1000 fiG) to about 188 MHz and 
680 MHz, respectively. Given the uncertainty of the mag- 
netic field strength during the synchrotron loss phase, the 
interstellar magnetic field in the area of S147, as well as 
the amount of shock compression (with or without cool- 
ing effects), this model cannot be ruled out as a possible 
explanation of the observed spectral bend. 

We conclude, in principle, both the compressed Galactic 
magnetic field scenario and early high synchrotron losses 
with subsequent adiabatic expansion could explain the ob- 
served turnover. The observed difference of Aa ~ 1 between 
the high-frequency and the low- freque n cy sp ectral indices 
is larger than the value given bv iTand (|1984D . but smaller 
than that from synchrotron losses with no further injec- 
tion of relativistic electrons after the first 100 years (see 
Fig. 2.12 in .Pacholczvk ..19771 

5. The polarization towards S147 

A6 cm and All cm polarization maps are shown in the lower 
panels of Figs. [1] and [2l The polarization features seen in 
the two polarization maps are complex, but very similar at 
both wavelengths. They have no counterpart in the total 
intensity map except for some sections of the S147 shell. 
Obviously, the interstellar diffuse emission component is 
responsible for most of the observed polarization features. 
The similarity of polarized structures at both wavelengths 
suggests that Faraday rotation along the line of sight is very 
small. 

We extracted the polarization intensity map of the same 
region (Fig. [TT|) from the A21 cm "EfF elsberg Medium 
Galactic Latitude Survey" (EMLS) (Re ich et all I2004D . 
This map also includes the large-scale polarized emis- 
sion component at an absolute zero level from the po- 
larization survey mad e with the DRAG 26-m telescope 
(jWolleben et al.( l2006( l. The polarized emission along the 
bright shell of S147 is very low. The shell of S147 might 
contribute polarized emission, which adds to the Galactic 
emission forming minimum polarization. As an alternative 
S147 might depolarize the emission from larger distances. 

Large-scale emission component is missing in the po- 
larization data at A6 cm and All cm, which is known 
to severely limit any interpretation of polarization struc- 
tures cau sed by Fara day rotatio n effects in the iii terstellar 
medium (|Reichll2006( ). Following [SuneLalJ (|2007l ) we used 
the WMAP polarization data at 22.8 GHz to recover the 
missing large-scale structures at All cm and A6 cm. 

For this purpose, we convolved the WMAP U and Q 
maps to 49' and scaled to All cm and A6 cm with a tem- 
perature spectral index of /? = —2.8, equal to the tem- 
perature spe ctral index toward the Galactic anti-center 
(j Reich fc Rei ch 1988). We obtained a spectral index of the 
polarized emission toward S147 oi (3 = —3.0 ± 0.4 from a 
T-T plot between the A21 cm map and the WMAP 22.8 



Effcsbcrg X21 cti PI 




1E2° lai' 130'' 17S° '78 = 

Galaclic Long'tJde 



Fig. 11. A21 cm polarization intensity map of S147 from 
the Effelsberg 100-m telescope at an angular resolution of 
9.'4 (grey scale). The r.m.s. -noise is 8 mK Tb. Contours 
show total intensities (the strong point-like sources were 
subtracted) starting from 150 mK Tb and increasing in 
steps of 100 mK Tb. The bars show the orientations 
of (E-|-90°). The length of the bars is proportional to the 
polarized intensity. 

GHz map, which are both on an absolute scale. Although 
the uncertainty is large, the spectrum of the polarized emis- 
sion is consistent with that for total intensity, implying that 
the polarized emission at longer wavelengths is not much 
reduced by depolarization and that a low Faraday rotation 
occurs in the direction of the Galactic anti-center. This jus- 
tifies the usage of the better known temperature spectral 
index /3 = -2.8. 

The difference of the scaled WMAP data and the All 
cm and A6 cm maps, which were convolved to the same 
beam, are then added to the original U and Q maps at A6 
cm and All cm wavelength. This procedure needs, how- 
ever, to be modified in case of significant Faraday rotation 
along the line of sight in particular for the All cm map. 
We have also used the A21 cm (Fig. [TT]) convolved to a 
49' beam and the WMAP 22.8 GHz map to obtain the av- 
erage Faraday rotation over the entire area. We obtain a 
value of ~ -1-10 rad m~^. A low Faraday rotation is also 
supported by the pulsar PSR J0538-I-2817, associated with 
S147, which has a Faraday rotation of —7 ± 12 rad m^'^ 
(|Mitra et al.ll2003[ ). We conclude that the correction of the 
polarization angle at All cm is smaller than a few degrees. 
In the following, we assume that no modification of the 
straightforward procedure is necessary. 

As expected for a magnetic field orientation along the 
Galactic plane a large offset is found for Stokes Q (average 
^ 4 mK Tb at A6 cm), while the offset for Stokes U is small 
(average ~ -0.2 mK Tb at A6 cm). The polarized intensity 
maps at All cm and A6 cm, after restoration for an abso- 
lute zero level, are shown in Figs. [T^ and [T51 The strongest 
polarized emission is located outside the boundary of S147 
in both maps. The distribution of polarized intensity in 
the A6 cm map does not allow the identification of S147 
in general. Only a few filaments, e.g., l,b — 181?6,— 1?0 
and l,b = 181?2, — 2?0, show an association between total 
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Fig. 12. Polarized All cm emission in the area of S147 
after restoration for an absolute zero-level using WMAP 
22.8 GHz data (grey scale). The angular resolution is 4.'4. 
Contours show total intensities from the Effelsberg All cm 
map (with strong point-like sources subtracted). The bars 
show the orientation of the magnetic field (E-|-90°) in case 
of ignorable Faraday rotation. The length of the bars is 
proportional to the polarized intensity. 
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Fig. 13. Same at Fig. [12] but for Urumqi A6 cm map. The 
angular resolution is 9.'5. 

intensity filaments and depressions in polarized intensity. 
In the All cm map, several filaments are clearly visible as 
depressions in polarized intensity. At both frequencies, the 
appearance of S147 against the large-scale polarized emis- 
sion component of the Galactic background is similar to its 
appearance at A21 cm (see Fig. [TT]), but the "depolariza- 
tion" effect is less serious at higher frequencies. 

The interpretation of the restored polarization data re- 
quires modeling of the diffuse Galactic polarized emission 
component including Faraday rotation effects (Reich 2006). 
However, the strong eastern filamentary structure is a re- 
gion where polarized emission from the SNR adds to the 



Galactic emission. It produces a depression in the polar- 
ized intensity map because the magnetic field along the 
filaments is oriented almost perpendicular to the general 
Galactic magnetic field direction. The strong polarized fil- 
ament structure is visible in the original polarization map 
at a relative zero level (Fig. [1] and Fig. [2]). 

5.1. Rotation measures of the eastern filaments 

For the calculation of rotation measures (RMs) of S147 we 
have to remove the diffuse Galactic polarized emission com- 
ponent. We noticed that the polarized emission from the 
eastern filament is strong and outstanding against the back- 
ground in the original polarization map. We calculated the 
RM map from the polarization angles of the original All 
cm and A6 cm polarization maps after both maps have been 
convolved to a common angular resolution of 10'. Fig. [HI 
shows the RM distribution (n=0) for the eastern filamen- 
tary region. We have excluded pixels with a polarized in- 
tensity below 1.40 mK T^ at All cm or 0.75 mK T^ at A6 
cm. 
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Fig. 14. The RM distribution for the eastern filamen- 
tary emission region of S147 derived from the original All 
cm and A6 cm polarization maps. Contours show polar- 
ized intensities from the Urumqi A6 cm map starting at 
2.5 mK Tb and running in steps of 0.8 mK T^. The posi- 
tion of the extrag alactic sour ce 0539-I-266B with a RM of 
-75 ± 4 rad m'^ (|Kimlll988f) is indicated. 

In principle, we require polarization maps at a minimum 
of three frequencies to determine the unambiguous RMs. 
However, we cannot use the A21 cm map for this purpose 
because the large-scale diffuse structure dominates and the 
filaments are not clearly visible. RMs calculated at A6 cm 
and All cm have an ambiguity of ±n x 362 rad m^^ (n=0, 
±1, ±2, ...). Any value of n other than zero seems unlikely 
toward the Galactic anti-center because the direction of the 
local magnetic field is almost perpendicular to the line of 
sight and thus small RM values are expected. 



10 



L. Xiao et al.: Spectral properties and the magnetic field of the SNR S147 



The average RM in the east filamentary region is found 
to be —70 ±8 rad m~^, which is almost identical to the RM 
=—75 ±4 rad of an ext ragalactic s ource, 0539+266B, 
shining through this region pCMll988D . Outside the fila- 
mentary shell, the RM excess caused by S147 is likely rather 
small, as indicated by the smooth distribution of polarized 
intensity and polarization angles in Fig. 1121 As was men- 
tioned previously the average Faraday rotation over the en- 
tire area is low. Therefore, the intrinsic magnetic field vec- 
tors of the eastern filaments are very close to the magnetic 
field direction shown in Fig. [TJ 

5.2. Magnetic fields of the eastern filament 

We estimate the magnetic field strength in the east- 
ern filamentary region by assuming energy equipartition 
between the magne tic field and electrons and protons 
(|Furst fc Reichll2004D : 

i?,ni„ « 10 • $ • R{deg) -6/7 . dikpc) -^/^ . S.GH.iJy) (1) 

We integrated the filamentary region and obtained a flux 
density at A6 cm of 5.1 Jy. The flux density at 1 GHz is 
extrapolated using the spectral index of a = —0.35. The 
radius of the radiating region is assumed to be 0?5 and the 
fraction of the radiation volume $ is ~ 1. For distances 
of 0.8 kpc, we obtain an estimate for the magnetic field 
strength of about 36/iG. Using the RM value —70 rad m^^ 
for the eastern filamentary region, we estimated the line of 
sight component of the intrinsic magnetic field. Assuming 
an electron density of 100 cm^-^ wit hin the filaments and 
their thickness to be 1" 0.07 pc (jFesen et al.lll985l ). we 
calculate « 35x RM/ng ^ 24/iG. Both values are of the 
same order. Magnetic fields in these filamentary structures 
are typically 30^G. 

6. Summary 

We present new sensitive maps of S147 for the total and po- 
larized intensity at A6 cm and All cm. The A6 cm maps are 
the first complete ones obtained so far. Our new measure- 
ments confirm the turnover of the integrated spectrum of 
S147 at a frequency of about 1.5 GHz. The measured struc- 
tures have been decomposed into filamentary structures 
(smaller than 18') and large-scale diffuse emission compo- 
nent. The T-T plots between 2639 MHz and 4800 MHz 
for these two emission components prove that the turnover 
is attributed entirely to the diffuse emission component, 
which has a spectral index of a ^--^-1.35. The filamentary 
emission component has a spectral index of a ~-0.35 and 
can be traced up to 40.7 GHz on WMAP maps. 

We discussed several possible mechanisms to explain the 
turnover of the diffuse emission component. We can ex- 
clude the scenarios of two populations of electrons as well as 
the diffuse shock accelerated of electrons. However, a com- 
pressed magnetic field is a possible explanation. The SNR 
shock wave compresses the local magnetic field and shifts 
the turnover in the Galactic radio spectrum from about 
400 MHz to about 1.5 GHz. High synchrotron losses, dur- 
ing the early phase of the SNR, would cause a bend at a 
rather high frequency, but subsequent expansion of S147 
would shift the frequency toward about 1.5 GHz as well. 

S147 is a faint source located on the near side of the 
Perseus arm. Polarized emission from S147 is confused with 



large-scale Galactic emission, except for the bright filamen- 
tary regions. The RM in the eastern filamentary region is 
about —70 rad m^^. The magnetic field in the filamentary 
region has a strength of typically ~ 30/iG. 
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